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ABSTRACT 

We present measurements of dust reddening using the colors of stars with spectra in the Sloan Digital 
Sky Survey. We measure reddening as the difference between the measured and predicte d colors of 
a star , as derived from stellar parameters from the SEGUE Stellar Parameter Pipeline ()Lee et al\ 
l2008al) . We achieve uncertainties of 56, 34, 25, and 29 mmag in the colors u — g, g — r, r — i, and i — z, 
per star, though the uncertainty varies depending on the stellar type and the magnitude of the star. 
The spectrum-based reddening measurements confirm our earlier "blue tip" reddening measurements 
jSchlafly et al. 2010, SIO), finding reddening coefficients different by -3%, 1%, 1%, and 2% m u — g, 
g — r, r — i, and i — z from those found by th e blue tip m e thod, after removing a 4% nor malization 
differe nce . These results pref er an Ry = 3.1 lFitzpatri"c3 ()1999l F99) reddening law to lO'Donnelll 
(Il994f) orlCardelli et al\ (|1989| ) reddening laws. We provide a table of conversion coefficients from the 
ISchleeel et al\ (|1998l . SFD) maps E{B — V) to extinction in 88 bandpasses for 4 values of Ry, using 
this reddening law and the 14% recalibration of SFD first reported by SIO and confirmed in this work. 

Subject headings: dust, extinction — Galaxy: stellar content — ISM: clouds 



1. INTRODUCTION 

Dust is composed of heavy elements produced by the 
nuclear burning of stars. These heavy elements are blown 
out of the stars in winds and explosions, and are re- 
processed i n the interste llar medium to eventually form 
dust grains (|Drainell2003( ) . The dust scatters and absorbs 
light, especially in the ultraviolet through infrared, ac- 
cording to the dust reddening law. The dust also emits 
photons thermally in the far-infrared. Accordingly, map- 
ping the dust is a central problem in astronomy. 

In previous work w ith the blue tip of the stellar locus 
(|Schlaflv et al\ [2010l . SIO), we examined the reddening 
law and the accuracy of the 'Sc hlegel et oil (IT998I . SFD) 
dust map usi ng photom etry from the Sloan Digital Sky 
Survey ( York et al.\ [200 0. SDSS) and the uniformity of 
the color of the blue tip of the stellar locus over the sky. 
This blue tip work recommended a 14% recalibration of 
the SFD dust map in the sense that E(B -V) = 0-86 ■ 
E{B — F)sFD, and a preference for a iFitzpatrickl (|1999L 
F99) reddening law over other reddening laws. In this 
work we set out to test that result using an independent 
set of data. 

The SDSS stellar spectra provide an independent test 
of reddening. Stellar spectra sensitively test reddening 
because the broadband photometry of a star is almost 
completely determined by three parameters: tempera- 
ture, metallicity, and gravity. These paramaters can be 
determined using only line information in the spectra, 
allowing the intrinsic broadband colors of the star to be 
predicted independently from the observed colors of the 
star. Dust intervening between us and the star, however, 
will shift the observed colors relative to the intrinsic col- 
ors. The difference between the predicted intrinsic col- 
ors and measured colors constitutes a measurement of 
the reddening to that star. This method is broadly sim- 
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ilar to that of iPeek and Graves! (|2010f ) and iJones et~al\ 
()2011|) , which use SDSS galaxy and M-dwarf spectra, re- 
spectively, to similar effect. 

The eighth data re lease of the SDSS has spectra for 
about 500,000 stars (lAihara et a l. 2011). The SEGUE 
Stellar Parameter Pipeline ^Lee'e t al. 2008a, SSPP) has 
uniformly processed these spectra to measure the tem- 
perature, metallicity, and gravity of each of these stars 
using a variety of methods, including ones that are inde- 
pendent of the observed colors of the star. We use one 
such method to predict the intrinsic colors of each star. 
The difference between the observed colors and intrinsic 
colors is used as a measurement of reddening to each star. 
These reddening measurements are then used to test the 
calibration of SFD and the reddening law. 

We test the calibration of SFD and the reddening law 
by comparing our reddening measurements with the pre- 
dictions from SFD and a reddening law. A reddening 
law predicts, in each color a — b, the reddening E(a — b) 
relative to the reddening in some reference color, usually 
B -V. The SFD dust map predicts E{B - V), effec- 
tively giving the normalization of the reddening law in 
any direction on the sky. We use our reddening measure- 
ments to find Ra-b = E{a — b)/E{B — F)sfd for each 
of the SDSS colors, testing both the reddening law and 
the SFD normalization. We also find the ratios of the 
Ra-b, which test the reddening law without using SFD 
as a reference. 

This work additionally lays the foundation for future 
investigations of the three-dimensional distribution of 
dust. Because we have stellar parameters for each of 
the stars, we can obtain accurate absolute magnitudes 
and distances. The SDSS targets both blue-horizontal- 
branch (BHB) stars and M-dwarf stars, in principle per- 
mitting the dust to be studied over a wide range of dis- 
tances. In this work, however, we have focused on the 
two-dimensional distribution of the dust. 

The paper is organized as follows: in §21 we present 
the data sets used in this work. In Sj3l we describe our 
method for transforming SSPP stellar parameters and 
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SDSS photometry into reddening measurements. In 21 
the reddening measurements are presented, cahbrated, 
and mapped. In ^ and Sj6l we discuss the imphcations 
of these results, especially for the reddening law, and 
conclude. 

2. DATA 

To measure the reddening to a star, we require its ob- 
served colors and a prediction for its intrinsic colors, as 
derived from stellar parameters for that star. The SDSS 
provides all of this information. In this section, we first 
describe the SDSS imaging, which provides the observed 
colors, and the SDSS spectroscopy and SSPP, which 
provide the stellar parameters. We then present the 
MAR CS grid of model atmospheres ()Gustafsson et al\ 
I2008D . which is used to connect the stellar parameters 
to colors. Finally, we describe the cuts we perform on 
the full set of SDSS stars with spectra to get the sample 
of stars used in this work. 

2.1. The SDSS 

The eighth data release of the SDSS provides uni- 
form, contiguous imaging of about one third of the sky , 
mostly at high Galactic latitudes (jAihara et a/.|[20Tll ). 
The SDSS imaging is performed nea rly simultaneously 
in five optical filters : u, g, r, i, and z (|Gunn et 
iFukugita et aLl[l996l ). The photometric pipeline has uni- 
formly reduced data for about 10^ stars. The SDSS is 
95% complete up to magnitudes 22.1, 22.4, 22.1, 21.2, 
and 20.3 in the bandpasses u, g, r, i, and z, which have 
central wavelengths of about 3600, 4700, 6200, 7500, and 
8900 A. We use SDSS data that have been photometri- 
cally cali brated according to the "ub ercalibration" pro- 
cedure of lPadmanabhan et al\ ()2008l ). 

The technique presented here is sensitive to the SDSS 
photometric calibration, and so we briefly describe the 
basic units in which the SDSS data is calibrated: the 
run and camera column ("camcol"). The SDSS takes 
imaging data in runs that usually last for several hours. 
These runs observe regions of sky that are long in right 
ascension and narrow in declination. Each run is com- 
posed of observations from six camcols: the six sets of 
five CCDs (one for each of the photometric bands) that 
fill the SDSS focal plane. 

The SDSS photometry is used to select targets for 
follow-up SDSS spectroscopy. Spectroscopy is performed 
with two multiobject double spectrographs, which are 
fed by 640 fibers. The spectra cover a wavelength range 
of 3800 - 9200 A at a resolution of A/AA ~ 2000. The 
stars used in this work were targeted for spectroscopy 
either as part of the main SDSS survey or as part of 
the Sloan Extension for Galactic Understanding and Ex- 
ploration (Ya nnv et al\ 12009 '. SEGUE) and its succes- 
sor, SEGUE-2. Targets from the main survey are ob- 
served until they reach a target signal to noise (S/N) 
of about (S/N)^ > 15/pix for stars with magnitudes 
g = 20.2, r = 20.25, i = 19.9. SEGUE targets are 
exposed long enough to achieve (S/N)^ > 100 at the 
same dept h, so that st ellar parameters can be measured 
([Abazaiia n a/.ll2009t ). 

The SDSS has a number of programs targeting different 
stellar types for spectra. Figure [T] shows a ugr color-color 
plot of the different types of stars u sed in this work, to - 
gether with the overall stellar locus. lYannv et al\ (|2009( ) 



gives a detailed description of the criteria that define 
these targets. The stars used in this work cover most of 
the stellar locus, except for the reddest stars, for which 
stellar parameters are poorly measured. Table [1] gives 
the names of each of the programs targeting stars used 
in this analysis, as well as the number of stars used. In 
addition to stars targeted for specific science goals, we 
also use the reddening and spectrophotometric standard 
star target types. These stars are F stars used for cali- 
bration purposes in the SDSS spectroscopic pipeline. 



f All SDSS siars 
f PHO 

RED 
f SEGUE1_BHB 

SEGUE1_FG 
^ SEGUEI^GD 
f SEGUE1_KD 
f SEGUE1_KG.' 




2.0 



0.0 



-0.5 



f All SDSS stars 
f SEGUE1_LM 

SEGUE1_MPMSTO 
f SEGUE2_LKG 

SEGUE2_LM 
f SEGUE2^MSTO 
f SEGUE2_PMKG 
f STAR_BHB •> . 




-0.5 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



3.0 



Figure 1. Colors of stars with SDSS spectra used in this work. 
The black contours give the number density of all point sources 
detected by the SDSS at greater than 20(t confidence in the ugri 
bands. The colored contours give the number densities of stars 
with SDSS spectra, targeted according to the program labeled in 
the legend. In each panel the horizontal axis gives u — g and the 
vertical axis gives g — r, in mags; the figure is split to make the 
contours more legible. 



2.2. The SSPP 

The Sloan Extension for Galactic Understanding and 
Exploration is a spectroscopic and photometric extension 
to the SDSS designed to provide insight into the structure 
and history of the Galaxy (Yannv et al. 2009) . As part of 
this program, the SSPP was developed to uniformly pro - 
cess the stellar spectra from the SDSS (ILee et aLII2008al) . 
This work uses the eighth data release of the SDSS, which 
includes an updated version of the SSPP, and spectra 
from SEGUE-2 in addition to SEGUE. The SSPP pro- 
vides multiple estimates of the spectral type, metallicity, 
gravity, and radial velocity for all of the stars with SDSS 
spectra. The accuracy of these parameters has been 
repeatedly tested by comparison with high resolution 
spectroscopy and by using rneasurements from globular 
and open clusters ('Lee et aLl |2008bl : IAllende Prieto efan 
[2008; .Smolinski et aL.2010, ). 
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Table 


1 




SDSS spectral targets 


-L ELrgct typo 


^ of stars 


Description 


PHO 


15894 


Spectrophotomctric standards 


RED 


14353 


Reddening standards 


SEGUE1.BHB 


12603 


BHB stars 


SEGUE1.FG 


5065 


EG stars 


SEGUEl.GD 


44495 


G dwarfs 


SEGUE1.KD 


13459 


K dwarfs 


SEGUE1.KG 


16012 


K giants 


SEGUE1.LM 


22273 


low nietallicity stars 


SEGUE1.LOW.KG 


2606 


low-latitude K giants 


SEGUE1.LOW.TO 


6432 


low-latitude turn-off stars 


SEGUE1.MPMSTO 


26885 


metal-poor F stars 


SEGUE2.LKG 


19634 


K giants 


SEGUE2.LM 


13785 


low nietallicity stars 


SEGUE2.MSTO 


31767 


main sequence turn-off stars 


SEGUE2.PMKG 


10566 


K giants 


STAR.BHB 


5667 


Main survey BHB stars 


Total 


261496 





Note. — The target types and the number of stars used for stars 
considered in this work. The low-latitude targets are selected according 
to targeting algorithms designed to s elect stars in highl y-reddened low 
Galactic latitude regions. SEGUE-1 I IYannv et a;j2009f) and SEGUE-2 
provide most of the targets for this work. The main survey, however, 
provides the STAR_BHB targets, as well as most of the spectrophomc- 
tric and reddening standards. 

The SSPP provides a number of different estimates for 
tlie temperature, metallicity, and gravity of each star, 
and a single composite estimate. We cannot use the com- 
posite estimates produced by the SSPP because some of 
the SSPP parameter estimates rely on the photometry of 
the stars, which is affected by the dust. Instead we rely 
only on methods that do not make use of any photometry. 
These are the N GSl, kil3, ANNRR, and ANNSR methods 
(|Lee et al\ 



2008a'). For this work we have used the ANNRR 



estimator ( Re Fiorentin et alll2007[ ). which uses a neural 



net trained on continuum-normalized SDSS spectra to 
make parameter estimates. However, the four methods 
produce compatible stellar parameters and the results of 
this work are insensitive to the method used. 

2.3. MARCS Model Atmospheres 

The stellar parameters derived from the SSPP are 
transformed into predicted ugriz colors using the 
MAR CS grid of model atmospheres ()Gustafsson et al\ 
120081) . Each model includes estimates of the surface 
flux at a resolution of A/AA of 20,000, and covers 
the range 1,300-200,000 A. The grid covers an exten- 
sive range of stellar parameters: 2500 K < Tcff < 
8000 K, -5 < [Fe/H] < 1, and -1 < log^ < 5.5, 
as well as additional parameter combinat ions involv- 
ing a-enhanccmcnt and microturbulence. lEdvardssonl 
([2008) and Plcz ( 2008) have verified the accuracy of syn- 
thetic broadband colors predicted from this grid. For 
log 5 < 3, we use the MARCS spherical stellar atmo- 
sphere models, while for log g > 3, we use the plane paral- 
lel models. We have expe rimented w ith using other grids: 
the iMunari et al\ ()2005D grid, the iCastelli and Kuruc3 
(12004^ grid and the NG Sl grid used internally in the 
SSPP. The IMunari et aU (|2005) grid predicts bluer col- 
ors than the other grids owing to the grid's not includ- 
ing "predicted lines," faint lines which h ave not been in- 
dividually detected ()Munari et a/.l [20051) . Nevertheless, 
because we calibrate the predicted colors to match the 
observed colors (|]4?2|), our final results are insensitive to 



the choice of grid, even when the raw predicted colors 
are substantially discrepant from the observed colors. 

2.4. Selection Cuts 

The SDSS targets a wide variety of stars for spec- 
troscopy, including exotic types for which the spectral 
parameters are difficult to measure. We are interested 
in the dust, and so want to select stars for which the 
stellar parameters are well understood and accurately 
measured. Accordingly, we exclude from this analysis 
stars targeted to have unusual colors and very red stars 
with dereddened g — r > 1 mag. Likewise, we exclude 
any white dwarf targets. Finally, we exclude any ob- 
jects with stellar parameters marked as unreliable by the 
SSPP (according the parameters' indicator variables). 

The remaining target types and number of stars used 
in each type are listed in Table [TJ We find it convenient 
to occasionally divide these target types into five classes. 
These classes are the standard stars, which were targeted 
as reddening standards or spectrophotomctric standards; 
the FG stars, targeted to have spectral types F or G; the 
BHB stars; the K stars, targeted to have spectral near 
K; and the other stars, which include the low-metallicity 
and low-latitude targets. 

We have chosen not to impose a cut on the S /N of the 
spectra. We have varied cuts on S/N from to 50; the 
resulting best fit reddening coefficients presented in this 
work vary by only 0.5% in that range. 

The total number of stars used following these cuts is 
261,496. 

3. METHODS 

We predict the colors of stars using the stellar param- 
eters from the SSPP and the MARCS grid. We compare 
these predicted colors to the measured colors from the 
SDSS imaging to derive the reddening to each star. 

The predicted broadband magnitudes nib are given by 
integrating each synthetic spectrum over th e SDSS sys- 
tem t hroughput for each of the five bands b (|Gunn et al\ 
|1998() , according to 



rrib 



-2.5 log 



JdXS{X)Wb{\) 
JdXZ{X)Wb{X) 



(1) 



Here the source spectrum S has units photons/s/A, the 
system throughput Wb is unitless, and the AB magnitude 
reference spectrum Z has the same units as S. The AB 
magnitude reference spec trum Z is simply a flat spec- 
trum with = 3631 Jy (|Oke and Gunn|[T983l ). 

The magnitude of a source in a single band predicted in 
this way is not useful because it depends on the unknown 
radius and distance of the source. However, in this work 
we consider only colors, for which the distance and ra- 
dius dependence cancels. We could in principle use the 
stellar parameters and stellar evolutionary tracks to de- 
termine the radius of the star, and so derive distances in 
conjunction with reddenings, but we defer this to future 
work. 

By computing the magnitudes corresponding to each 
synthetic spectrum, we construct a synthetic grid of mag- 
nitudes. The magnitudes for each star are then predicted 
by linearly interpolating off this grid. The differences be- 
tween the measured colors and predicted colors give red- 
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dening estimates. We define this reddening measurement 
in the color a — 6 to be Aa-t = ~ &)obs — (a — 5)pi.c. 

4. RESULTS 

The reddening estimates Aa„b are used to measure 
the dust reddening law and test dust properties. In this 
section, we describe the performance of these reddening 
measurements at high Galactic latitudes f ti4.ip . demon- 
strating the need for an empirical calibration to improve 
the measurements. We then perform that calibration 
( §4.2p . Finally, we present maps of the calibrated red- 
dening measurements ( §4.3p . 

4.1. Performance at High Galactic Latitudes 

We test the quality of the A by restricting to the re- 
gion of sky with |6| > 50° and E{B - V)sfd < 0.04. 
In this region of low extinction, excursions of A from 
zero should be due primarily to the statistical uncertain- 
ties in the photometry and stellar parameters. However, 
additional systematic effects can cause A ^ 0. These 
include any mismatch between the synthetic spectra and 
the real spectra and biases in the SSPP-derived stellar 
parameters. 

The performance of A for this sample of stars is given 
in Table [U for several classes of star. The rows marked 
"raw" give the means and standard deviations of the 
uncalibrated A relevant here. In all colors, the mean 
A is within about 30 mmag of 0, comparable with 
the 10-20 mmag unce rtainty in the SDSS zero points 
(jAbazajian ltd\\2(M t. In r — i and i — z, the scatter in 
A is 20-30 mmag, comparable to the SDSS photometric 
uncertainty. 

However, in the colors u ~ g and g — r, the scatter in 
A is somewhat larger than expected. The expected un- 
certainty in A„_g and Ag_r is about 40 and 30 mmag, 
respectively, including both the photometric uncertainty 
and the uncertainty in the stellar parameters. However, 
the scatter in A in this sample of stars is 92 and 51 mmag 
in the colors u — g and g — r. Figure [2] (left panels) shows 
that additionally there are trends in A with tempera- 
ture, metallicity, and gravity, especially in the color u—g. 
Because we expect that typical stellar metallicity is cor- 
related with the dust column to that star, these trends 
could bias this analysis and must be removed. Moreover, 
these trends contribute to the scatter in A„_g and Ag-r- 

We do not fully understand the source of these trends. 
The two most likely candidates are mismatches between 
the synthetic stellar spectra and real stellar spectra, and 
biases in the SSPP stellar parameters. We have tested 
for mismatch between the synthetic and observed spectra 
by comparing with other synthetic grids. Comparing the 

t'ds of [Gustafsson et al. (2008), Castclli and Kurucz 
)03), and ILee et ali (2008a) indicates that the grids 
can disagree with one another by about 50 mmag in u — g 
and g — r. The larger effect, however, is probably due to 
biases in the SSPP stellar parameters. The data release 
8 version of the SSPP introduced substantial changes to 
metallicity estimates fo r stars ar ound solar metallicity, by 
about 0.4 dex (Smohn ski et all 12010.). Colors predicted 
using parameters from the data release 7 version of the 
SSPP were too blue by as much as 200 mmag in it — 5. 

However, we calibrate the predicted colors to re- 
move the trends and to render our results insen- 
sitive to these sources of systematic errors. We 



have repeated the full analysis using four different 
grids (iGustafssoii et a/.l 120081: iCastelh and Kurilc3 120041: 
ILee et a/.ll2008al: iMunari et a/.ll2005f) . The final results 
of this work were unchanged at the 1% level. We have 
also repeated the analysis using the data release 7 and 8 
versions of the SSPP, which changed the final results by 
at most 2%. 

4.2. Calibration 

We remove these trends by calibrating the syn- 
thetic photometry so that A « on the high- 
latitude sample of stars where reddening is unimpor- 
tant. This procedure has b e en do ne before; for in- 
stance, iFitzpatrick and Massal ()2005[ ) also calibrate syn- 
thetic photometry to match observed photometry. In this 
work, we perform the calibration by fitting an empirical 
curve to A as a function of temperature, metallicity, and 
gravity, and adjusting the synthetic photometry to re- 
move this curve from A. 

To properly model the trends in A, we should find 
the function / of temperature T^s, metallicity log Z, and 
gravity log 5 that most accurately predicts the colors. 
To be more precise, / should be chosen to minimize the 
distance between A and / in color space, considering 
the covariance matrix for the observed colors and for the 
predicted colors. We however deemed this computation- 
ally intractable and unnecessary, given the wide range 
of temperatures (4000 K < Tcff < 8000 K) , metallicities 
(—3 < logZ < 0), and gravities (2 < log 5 < 5) avail- 
able and the small estimated uncertainties in the stel- 
lar parameters for the objects we look at (ctt < 150K, 
o'logZ ^0.1 dex, CTiogg ^0.1 dex). We therefore ignore 
the covariance between the stellar parameters and A, 
and rely on the fact that the stars cover a broad range 
of parameters (much broader than the uncertainties in 
those parameters) to guarantee that the bias introduced 
by this is small. 

To maximize the effective range of stellar parameters, 
we need to a set of stars for calibration that cover param- 
eter space evenly without overweighting stars of a par- 
ticular type. Accordingly we choose stars from the ugr 
color plane, capping the maximum number of stars cho- 
sen from any location in the plane. This guarantees that 
all of the available stars are used in areas of color space 
where few stars were targeted, but limits the number of 
stars used in areas where many stars were targeted. The 
primary consequence of this selection is to reduce the im- 
portance of the many SEGUE1_GD targets chosen from 
a narrow range in temperature. 

With this set of stars in hand, we find the function / 
that minimizes 

i ^ ' ^ 

for each color c, where i indexes over stars, A gives the 
measured color minus the predicted color, and x gives 
the stellar parameters for the stars. The uncertainty a 
is computed from the photometric uncertainty and the 
uncertainty in the predicted colors derived from the un- 
certainty in the stellar parameters. The values A here 
are corrected for reddening according to SIO, though as 
E{B — F)sFD is less than 0.04 for this sample this is of 
little importance. We use for / a fifth order polynomial. 
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Table 2 

Bias in Predicted Magnitudes 



Target class 



Standards (raw) 


22 


51 


3 


33 


28 


22 


17 


25 


Standards (cal) 


-3 


46 


-2 


31 


-1 


22 


-1 


24 


FG (raw) 





82 


-2 


42 


30 


28 


20 


35 


FG (cal) 


1 


72 


-2 


37 





27 





34 


BHB (raw) 


64 


83 


-1 


50 


36 


32 


33 


42 


BHB (cal) 


24 


90 


-13 


46 





31 


-2 


40 


K (raw) 


10 


114 


-24 


57 


29 


32 


14 


32 


K (cal) 


5 


73 


4 


35 





27 





28 


Other (raw) 


13 


88 


20 


47 


34 


28 


18 


31 


Other (cal) 


-16 


68 


8 


38 


3 


27 





30 


All (raw) 


12 


92 


-6 


51 


30 


29 


18 


33 


All (cal) 


1 


71 





37 





27 





31 



Note. — Th e me an (A) and standard deviation (c r) of A in mmag for different 
types of stars ( ^2.2| |. before and after calibration ( q4.2l l. in each of the SDSS 
colors. Rows marked "(raw)" or "(cal)" indicate whether the means and standard 
deviations are given before or after calibration, respectively. After calibration, the 
mean of A is close to zero, and the scatter in A is reduced, especially in the color 
u- g. 





Figure 2. Color difference A in mags with temperature Teg , metallicity log Z, and gravity log g for the colors u — g, g — r, r — i, and i — z. 
The left hand panels show A before calibration, while the right hand panels show A after calibration. Before calibration, the u — g color 
is worst predicted, and small but easily detectable trends with temperature exist for all of the colors. Note the different scale for u — g. 



though we have also used a second order polynomial with 
negligible effect on the final results of this work. 

Figure [5] (right panels) shows that the remaining trends 
in A are small. The rows in Table [2] marked "cal" show 
the means and standard deviations of A after calibration 
for several classes of star. The calibration brings the 
mean of A to zero. The calibration has the additional 
benefit of reducing the scatter in A to be much closer to 
the expected statistical uncertainty of about 40, 30, 25 
and 30 mmag in u — g, g — r — i, and i — z. 

To verify that the final results are insensitive to the 
details of the calibration, we have varied the order of the 
polynomial used and the sampling method used to select 
calibration objects. As long as the order of the polyno- 



mial is at least two, we achieve consistent results at the 
1% level. The sampling method is similarly unimportant 
except in extreme cases where the calibration stars do 
not cover a satisfactory range of stellar parameters. 

4.3. Reddening Map 

The calibrated measurements of A immediately per- 
mit the construction of a reddening map, which is shown 
in Figure [3] (left panels) . The right panels show the map 
after correction for dust according to SFD, using the pre- 
scription of SIO. The derived reddening map agrees well 
with that found in SIO, clearly identifying the same runs 
with bad zero points in u— g and the SFD-underpredicted 
region in the northwest of the north Galactic cap. In the 
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Galactic plane, the stars used for the reddening measure- 
ments may not be behind the entire dust column, and so 
it is unsurprising that their reddening is less than that 
predicted by SFD. Over most of the observed sky, how- 
ever, the residuals after correction for dust are within the 
uncertainties. 

As pointed out in SIO, we again find that the dered- 
dened south is somewhat redder than the dereddened 
north fi jS.Sp . Apart from this trend and the cloud in the 
northwest, the residual maps are reassuringly flat and 
reveal no trends in color with Galactic latitude or longi- 
tude. 

Figure [H shows t he residual g — r map of 
iPeek and Gra^ (|20lol ). This map is compared with our 
residual g — r map, smoothed to match the 4.5 degree 
resolution of the IPeek and Graves! ([20101 ) map. The two 
maps show several consistent large scale features. Both 
maps detect the SFD underprediction in the northwest, 
and additionally find that the southwest of the North 
Galactic Cap is somewhat redder than the southeast, by 
about 10 mmags. 

5. DISCUSSION 

The reddening measurements A permit tests of the 
reddening law and of the normalization of SFD. In this 
section we use A to measure reddening in the SDSS 
colors, simultaneously constraining the normalization of 
SFD and the reddening law, using two methods. In tj5.11 
we find the best fit values of Ra-b, the reddening per 
unit E{B - V)sFD in each of the SDSS colors. In ^Ol 
we find the best fit values of the ratios of the Ra-b, with- 
out using SFD. In §5.3| we discuss two possible sources 
of systematic error in the measurements: photometric 
bias from target selection and the bias introduced by in- 
cluding stars at different distances. In H5.4\ and we 
examine the variation of reddening in different regions of 
the sky and a curious difference in residuals between the 
Galactic north and south. 

5.1. Fitting the Reddening in the SDSS Colors 

We use the reddenings A to constrain Ra-b, the red- 
dening in the color a — b divided by E{B — V^)sfd- We 
measure the Ra-b by fitting A as a linear function of 
E{B — F)sFD, the slope of which gives Ra-b- We find 
that we can detect calibration errors in the SDSS using 
this technique. Accordingly, we perform a second similar 
linear fit, where additionally a constant zero point offset 
in each SDSS run is also fit. 

The Ra-b are given by the equation 

^a-b = Ra-bEiB ~ F)sFD + C (3) 

where C gives a constant offset which should give zero 
if the calibration of the reddening values in the high- 
Galactic-latitude, low-reddening is accurate when ex- 
tended over the entire footprint. We find Ra-b through 
least-squares minimization of 

^2 _ ^ ( ^a-b,i — Ra-bE{B — V)sFDA — 

where i indexes over stars, Ra-b and C are the fit pa- 
rameters, E{B — V)sFD,i is E{B — V)sFD in the direction 
of star i, and ai is the uncertainty in Aa-b.i, consider- 
ing both the uncertainty in the measured color and the 



Table 3 

B.a-b derived by fitting SFD 



Color 


Ra — b 




Ra-b (zp) 




SIO 


SFD 


u - g 


0.94 ± 0.02 


1.45 


0.98 ± 0.01 


1.41 


1.01 


1.36 


g - r 


0.98 ± 0.02 


1.19 


0.94 ± 0.01 


1.12 


1.00 


1.04 


r — i 


0.55 ± 0.01 


1.04 


0.56 ± 0.01 


1.01 


0.57 


0.66 


i — z 


0.44 ± 0.01 


1.02 


0.42 ± 0.01 


0.97 


0.45 


0.61 



Note. — Ra — b derived from fitting SFD. The column Ra — b gives 
the fit results without zero point offsets, while the column Ra-b (zp) 
gives the fit results with zero point offsets. The P^r degree of free- 
dom {x^ /N) is close to unity in the redder bands, and improves when 
zero points are fit (x^p/^)- The predictions for Ra — b according to 
ISchlafiv e7~an | |20ia SlO) and the combination of SFD and lO'DormelH 
are also given. The spectrum-based R^-h arc seen to closely 
match one another and the SIO values, while disagreeing with the pre- 
diction according to SFD and iO'DonneU & 19941) . 

uncertainty in the predicted color, propagated from the 
uncertainty in the stellar parameters. 

To render the fit insensitive to calibration errors in the 
SDSS, we also fit 

Aa-b = Ra-bEiB - V)SFD + C, (5) 

in an analogous least-squares sense, where r indexes over 
SDSS run number, which absorbs any zero point cali- 
bration errors in the SDSS runs. We call this the fit 
with zero point offsets. The results of this fit are sim- 
ilar to those given by Equation [31 Each fit is iterated 
and clipped at 3cr in each color, until finally all clipped 
points are removed and the fit is repeated with the same 
clipped set of stars in each color. 

The results from the fits are given in Table [3] The 
fits with and without zero point offsets agree to within 
4%. Both sets of fits give values consistently slightly 
smaller than the SIO values 4%). In SIO, we found 
that the best fit normalization of SFD varied over the 
sky by about 10%, and so this 4% normalization differ- 
ence may be due to the different sky footprint available 
to the SIO and spectrum-based analyses. However, when 
we tested this at high Galactic latitudes in the north 
where the two footprints overlapped, the results were in- 
conclusive because the uncertainties were similar in size 
to the 4% effect we observe. The per degree of free- 
dom (x^/dof) for the fits are quite good, especially in 
the redder colors; without including zero point offsets, 
we achieve 1.45, 1.19, 1.04, and 1.02 x^/dof in the colors 

~ 9, 9 ~ r — i, and i — z, respectively. 

These results confirm the preference for an Ry = 3.1 
F99 reddening law and a 14% recalibration of SFD, first 
reporte d in SIO. The pr edicted Ra-b according to SFD 
and an lO'Donnelll (Il994() reddening law are in disagree- 
ment with the measurements. 

Figure [5| shows the residuals to the fit without zero 
point offsets. The residuals are largely Gaussian, with 
a of 56, 34, 25, and 29 mmag in u — g, g — r, r ~ i, 
and i — z. This scatter is substantially better than the 
scatter we achieved in the calibration sample of stars 
(Table [2]) , because those stars were chosen to be widely 
distributed over the ugr color-color plane and therefore 
many of those stars had outlying photometry. 

The residuals as a function of E{B — 1^)sfd display a 
striking trend, shown in Figure [6] The residuals increase 
with E{B - V) for E{B - V) < 0.2, and then decrease 
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Figure 3. Map ot A in u — g, g — r, r — i, and i — z (left panels), giving the derived reddening to each star targeted by the SDSS. 
The right hand pane ls (marked "corr") give the same maps, after correction for the dust according to SFD and the coefficients Ra-b from 
ISchlafiv et ali 120101) . The color bars give the range of A in mmag corresponding to the grayscale figure. The horizontal axis is right 
ascension and the vertical axis is declination, both in degrees. 



thereafter. This result is consistent with the SIO resuh 
that the best fit normalization of the dust was larger 
for E{B -V) < 0.2 than it was for E{B - V) > 0.2 
by about 15%, which is about the size of the trend seen 
here. That said, the entire E{B — V) > 0.3 footprint 
for this analysis occurs in a relatively small area near 
the Galactic anticenter and a few SEGUE low latitude 
fields, and may not apply more generally. 

The formal statistical uncertainties in the fit results 
are small (~ 2 parts in a thousand). These uncertain- 
ties dramatically underestimate the true uncertainty in 
the fit, stemming from biases in the photometry ( §5.3.1|) . 
unmodeled calibration errors in the SDSS, and problems 
with using SFD as a template (SIO). The reported un- 
certainties in Table [3] are derived by a Monte Carlo sim- 
ulation with simple models for the calibration errors and 
SFD template errors. Calibration errors in the SDSS 
are modeled as normally distributed run zero point er- 
rors, with standard deviations of 20 mmags in u and 
10 mmags in g, r, i, and z . according to the es t imate d 
calibration uncertainties of iPadmanabhan et al\ ()2008[ ). 



We additionally allow the true zero point to change by 
at over the course of the night, where d gives the rate 
of change of zero point and t is the time of night. We 
choose a to be normally distributed with standard devi- 
ations of 20 mmags per 6 hours in u and 10 mmags per 
6 hours in g, r, i, and z. For the SFD template errors, 
we multiply SFD by a normalization which varies over 
the sky on 1 degree scales, with a mean of 1 and stan- 
dard deviation of 0.1. This is very roughly intended to 
simulate the uncertainty in the SFD dust temperature 
correction, but does not include the covariance between 
temperature and extinction that may be important. 

The resulting uncertainties are dominated by uncer- 
tainties coming from the calibration effects, with the zero 
point and change in zero point with time effects con- 
tributing similarly. 

5.2. Fitting Ratios oj Ra-b using A Color-Color 
Diagrams 

We can free the analysis from a dependence on SFD 
and its normalization by looking only at ratios of Ra-b- If 
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Figure 4. Color ma ps of the North Galact ic Cap vsi g — t 
after dereddening, from IPeek and Graved 1120101 ) and the stellar 
spectrum-based measurements of this work, smoothed to 4.5°. 
Some of the same structures are visible. The SFD underpredicted 
region in the northwest is obvious, but likewise in both sets of red- 
dening measurements the southwest and northeast are redder than 
the southeast. Both maps are scaled so that black is —10 mmag 
and white is 10 mmag. Axes give coordinates in right ascension 
and declination, in degrees. 
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Figure 5. The distribution of residuals x (Equation |4)l in mmag 
for the fit without zero point offsets. The residuals are largely 
Gaussian, though the wings are non-Gaussian and are probably 
caused by errors in SFD and the changing best-fit dust normaliza- 
tion. A Gaussian fit is overplotted, with standard deviation given 
by a in each panel. 

A has been appropriately calibrated, then in the absence 
of dust A should be consistent with zero. The presence of 
dust will spread A along a line in color-color space. The 



Table 4 

Color-color Fit Results 



Color 


Color-color fit 


SIO 


O'Donnell 


Rg — rj Ru — g 


0.99 ±0.02 


1.00 


0.77 


Rr — i/ Rg — r 


0.59 ±0.02 


0.57 


0.64 


Ri — zl Rr — i 


0.77 ±0.02 


0.79 


0.91 


Note. — 


Ratios of Ra—b derived 


from fitting 


A color-color 



diagrams. The spect rum-based Ra-b are seen to closely match the 
ISchlafly et ali II2010I . SIO) measurement, while disagreeing with the 
O'Donnell prediction. 

slope of this line in a pair of colors a — b and c — d gives 
the ratio Ra-b/ Rc-d and is determined by the reddening 
law. 

We find the best-fit the slope of the line in A in color 
space in a least-squares sense, considering the covariance 
between the Aa-b induced by the photometry and the 
color predictions from the stellar parameters. Figure [7] 
shows the result of this fit, together with th e predicted 
line from SIO and from an lO'Donnel l (1994) reddening 
law. Table |4] tabulates the best fit parameters for the 
ratios of Ra-b- The fits do a good job at r eproducing 
the SI O ratios for Ra-b, but disagree with the lO'Donnelll 
(|1994f) reddening law predictions. We achieve x^/dof of 
1.4 on average over all the stars we fit. On the more 
limited subset of standard stars we get x^/dof of 1.1. 

An important caveat applies to this analysis. The 
Aa-b for a star are correlated. We track the correla- 
tion induced by the uncertainty in the photometry and 
the uncertainty in the stellar parameters. However, any 
unknown effect that changes the predicted or measured 
colors of stars will additionally spread out A in color- 
color space, and there is a danger of interpreting that 
signal as signal from the dust. When instead we fit A 
as a function of E{B — T^)sfd, we have the advantage 
of using a template that is uncorrelated with the uncer- 
tainties in A and that is known to be tightly correlated 
with the dust. 

To verify the robustness of these results we have sepa- 
rately fit the standard stars, the F and G stars, the BHB 
stars, the K stars, and all of the target types we look at. 
The results were consistent at the 3% level. 

5.3. Sources of Systematic Error 

5.3.1. The Photometric Bias 

Most SDSS stars are targeted for spectra on the ba- 
sis of their photometry. Usually stars are targeted for 
spectra if they fall within a specific box in color space. 
Photometric uncertainty may scatter stars into the color 
box preferentially from one side if the box is placed in a 
region of color space where the gradient in stellar num- 
ber density with color is different on one side of the box 
than the other. This effect will lead to a bias in the 
photometry. 

Were this bias a constant for the entire survey it would 
be taken out by the calibration of A. However, most of 
the SDSS and SEGUE color boxes select stars on the 
basis of their SFD-dereddened colors. Thus, errors in 
SFD and its normalization can essentially shift the color 
box relative to the stellar locus and introduce a bias that 
changes systematically with E{B — U)sfd- 

We have analyzed the magnitude of the bias introduced 
by this effect by simulating the photometric bias for each 
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Figure 6. Residuals x (Equation [4J in mmag for the fit without run offsets, as a function of E{B — V) in mag. The four columns 
give the colors u — g, g — r, r — i, and i — z. The rows show the residuals for different types of stars: the standard stars, FG-type stars, 
BHB-type stars, K-type stars, other stars, and all of the stars used in this work. There is a consistent trend in the residuals to increase 
when E(B — V) < 0.2, and to decrease thereafter. 



of the SDSS and SEGUE target types. We take the SDSS 
stellar locus at high Galactic latitudes as a proxy for the 
true stellar locus. Errors are artificially added to the 
stars' magnitudes consistent with the SDSS estimates for 
the uncertainties to create a new observed stellar locus. 
The mean error of a star making it into a color box gives 
the photometric bias. The change in the photometric 
bias with E{B—V)sfd can be investigated by shifting the 
colors of the stars according to the error in E{B — V^)sfd- 
For the target types considered in this work, the pho- 
tometric bias is less than 0.02 mags, though for very red 
or unusual types it can be much larger. The potential 
bias introduced in the derived Ra-b is 2%. 

5.3.2. The Three Dimensional Footprint 

The stars targeted by the SDSS for spectra probe a 
range of distances. For the dwarf stars in the survey, we 
use the iJuric et ah ()2008i) "bright" photometric parallax 



relation to get approximate distances to stars from their 
colors. We find that the stars we analyze have distances 
from 1-6 kpc, with a long tail to greater distances. 

The SFD dust map was calibrated to match the red- 
dening of galaxies, which will be behind any dust in the 
Galaxy. If substantial dust exists beyond the distances 
probed by the stars in our data, the Ra-b found in this 
work will be less than the true values. However, we find 
compatible Ra~b for the K stars 1 kpc away) as for 
the standard stars in the sample 4 kpc away), indicat- 
ing that the dust is confined to within 1 kpc over most of 
the SDSS footprint. In future work we intend to study 
the three-dimensional distribution of the dust by prob- 
ing nearer distances with redder stars and by focusing on 
regions of lower Galactic latitude. 



5.4. Fits to Different Sky Regions 
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Figure 7. Color-color plots of A in mag. The data should have A consistent with z ero in the absence of dust. Reddening from dust 
spreads A along a line given by the reddening law. The predictions for this line for the lSchlafiv et all l|2010l . SIO) measurements and for 
an O'DonncU (1994) reddening law are given by the blue and green lines. The red line is the best fit to the data, and is consistent with 
the SIO measurements but not with the lO'Donnelll 1)19941 ) prediction. The oval in the lower right of the panels gives a typical uncertainty 
covariance ellipse for a point on the diagram. 



The A are sufficiently numerous and cover a sufB- 
ciently large region of the sky that they can be usefully 
cut into different subsets to study the variation in Ra-b 
over the sky. Figure [8] shows the results of such an analy- 
sis. The panels show Ra~b for different ranges in Galactic 
latitude, longitude, E{B — V)sfd^ and dust temperature 
according to SFD. The rightmost panel shows Ra-b as 
derived for stars observed with each of the SDSS cam- 
cols. All Ra-b are fit with zero point offsets to minimize 
sensitivity to calibration errors in the SDSS. 

The fits results show smaller scatter between the vari- 
ous sky regions. We recover the SIO result that stars with 
E{B — y)sFD < 0.1 prefer a larger normalization than 
other stars by ~ 10%; likewise in the north relative to 
the south, and possibly for the same reason, as the north 
has less dust than the south. Fits to the dust in different 
ranges of Galactic longitude and at different tempera- 
tures according to SFD are reasonably consistent. 

We can get some idea as to whether or not the filter 
differences between the SDSS camcols are important to 
this analysis by splitting the survey by camera column. 
We find very close (2%) agreement between the various 
camcols. This also serves to place an upper bound on 
the statistical uncertainty in the measurements. 

These results on the universality of the Ra-b are con- 
sistent with those found in SIO. In that work, we did a 
more detailed analysis of the variation in Ra-b over the 
sky and found the dominant mode of variation in Ra-b 
took the form of an overall 10% variation in the overall 
normalization of the Ra-b- We have not repeated that 
analysis here, but we try to account for the effects of a 
10% normalization variation in our uncertainties f iil5.ll) . 

5.5. North/South Color Asymmetry 

In SIO, we noted an asymmetry between the north and 
south in the dereddened color of the main-sequence turn- 
off, which we called the blue tip color. We were unable to 
tell whether this color difference was a calibration effect 
or a result of the south having an intrinsically different 
blue tip color than the north. Because the spectrum- 
based reddening measurements are sensitive to calibra- 
tion problems and not to the intrinsic color of the blue 



Table 5 

Color difference between North and South 



Method 


u- g 


g-r 


r — i 


i — z 


Spectrum 


2.3 ±7.2 


8.8 ± 1.5 


3.4 ± 1.9 


9.3 ± 1.6 


Blue tip 


7.6 


21.8 


7.2 


12.4 



Note. — Median color difference in mmag between the north 
and south, for 40 < |fe| < 70, using the spectrum-based and blue 
tip methods for measuring colors. A calibration or dust effect will 
create a color difference in both sets of measurements, while the 
presence of a stellar structure in the south would create a color dif- 
ference only in the blue tip measurements. The blue tip offsets in 
r — i and i — z are similar to those in the spectrum-based measure- 
ments, suggesting that these color differences do not come from 
a stellar structure. However, in g — r, the color difference could 
indicate different stellar populations in the north and south. In 
u — g, systematic uncertainty in the offsets prevent determination 
of whether the offset comes from calibration or is astronomical. 

tip, we can use them to break this degeneracy. 

Table [S] shows the differences in A between the Galac- 
tic north and south, compared with the blue tip color 
differences found in SIO. We find the offset in A for each 
target class, and report the mean of these offsets in Ta- 
ble [5] For the uncertainty in the means, we report the 
standard deviation of the offsets. The color differences 
are comparable in r — i and i — z, suggesting that calibra- 
tion errors in the SDSS are the source of the discrepancy 
in these colors. In g — r, however, the blue tip color dif- 
ferences are larger than the spectrum-based differences. 
In It — g, the uncertainty is too large to make any defini- 
tive statement. The difference in these colors could then 
be caused by a fundamentally different blue tip color in 
the north than in the south. This result is heartening 
because the color of the blue tip varies due to age and 
metallicity changes in stellar populations, which should 
primarily affect the u — g and g — r colors while leaving 
the r — i and i — z colors unchanged. 

6. CONCLUSION 

The SDSS provides a wealth of information for test- 
ing reddening. Using the SSPP stellar parameters, our 
spectrum-based technique achieves reddening measure- 
ments with empirical uncertainties of 56, 34, 25, and 29 
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Figure 8. Ra—b for various subsets of the full data for cacii of the SDSS colors, using fits with run offsets. Shown are Ra—b for the north 
{b > 25°), south {b < —25°), and plane (|fe| < 25°); Ra~b for four ranges of Galactic longitude; Ra-b for regions of increasing E(B — V)sfd; 
Ra-b for regions of different temperature; and finally Ra—b for the 6 SDSS camera columns. Overplotted in blue are the 'Schlafiy et ali 
(|2piO) values divided by 1.04 ( i|5.1l l: the thick gray line gives the fit over all the data. Error bars are slightly offset horizontally for legibility. 
The fitted Ra—b a-rc reasonably uniform over the sky, except most noticeably for a 10% normalization difference between the fit results in 
the Galactic north and south. 



mmag in the colors u — g, g — r, r — i, and i — z, compa- 
rable with the expected uncertainties from the photom- 
etry and stellar parameters. Spectra of a few blue stars 
on a sight line are sufficient to tightly constrain redden- 
ing provided that well-calibrated photometry is available 
and that enough spectra are available to calibrate the 
predicted colors to the measured colors. 

The use of individual stars with spectra additionally 
allows extinction to be studied as a function of distance. 
Stellar spectra permit good estimates of the intrinsic 
luminosity of a star in addition to its reddening, mak- 
ing three-dimensional studies of the dust feasible. The 
wide range of SEGUE target types — M and K dwarfs 
through BHB stars — will likely make such analyses espe- 
cially fruitful, at least at low Galactic latitudes where all 
of the stars are not behind the entire dust column. At 
high Galactic latitudes, we could in principle use large 
numbers of stars to test for reddening through interme- 
diate and high-velocity clouds. We defer to later work 
the attempt to extend this analysis into the third dimen- 
sion, having established the feasibility and accuracy of 
the method in two dimensions. The combination of large 
photometr ic and spectroscopic surv eys like Pan-ST ARRS 
i|Kaiser et at 2002) and LAMOST (|Su et aLlll998l ) seems 
particularly promising. 

These tests should also provide useful feedback to mod- 
els of synthetic spectra and to stellar parameter esti- 
mates. We look forward to incorporating the DR8 ver- 
sion of the SSPP into this analysis, and potentially pair- 
ing it with alternative synthetic spectral grids. However, 
we note that while this analysis provides an effective test 
of the colors predicted by synthetic spectral grids, never- 
theless it does not depend on the accuracy of these colors. 
The wide range of spectral types targeted by the SDSS 
permits the predicted colors to be well calibrated to the 
observed colors, rendering the final results independent 
of the color accuracy of the synthetic spectra. 

The spectrum-based reddening measurements give 
best-fit Ra^b that closely agree with the SIO values. Ac- 



cordingly, we have gained confidence in the F99 redden- 
ing law with Rv = 3.1 and normalization N — 0.78 
proposed in SIO. The variation in the best-fit normaliza- 
tion of the reddening law seen in that analysis remains 
a problem. The spectrum-based reddening tests permit 
the variation in best-fit normalization with extinction to 
be seen more clearly; when E{B — V)sf'd ^ 0.2, the 
best-fit normalization is about 15% higher than when 
E{B — V)> 0.2, though this conclusion relies on the rel- 
atively small fraction of the sky where E{B — V) > 0.3 
and SDSS data is available. 

Nevertheless, the agreement of the SIO and spectrum- 
based reddening measurements demands that the F99 
reddening law be used to predic t reddening over 
ICardelli et aJi (fT989t ) or lO'Donnelll ([1991 reddening 
laws. Appendix [S] gives the extinction per unit E{B — 
^)sFD predicted by these sets of measurements in 88 
bandpasses for 4 values of Ry- We have focused on high 
Galactic latitudes where we observe Ry = 3.1; other val- 
ues of Rv are provided only for convenience. Because we 
have seen that the best fit reddening law normalization 
varies over the sky and as a function of E{B — V)sf'd, 
it is possible that outside the SDSS footprint a different 
normalization might be preferable. However, the shape 
of the reddening law seems constant over the SDSS foot- 
print, and the normalization we suggest is unambigu- 
ously the best choice over the large area of sky covered 
by the SDSS. Therefore we propose that this reddening 
law and normalization become the default choice to be 
used in the absence of other information. 

Multiple sensitive, mutually consistent measurements 
of reddening over the SDSS footprint are now available. 
Extension of this work to larger areas of sky seems read- 
ily possible as Pan-STARRS and LAMOST data become 
available. These measurements will permit a next gener- 
ation dust map to be constructed and tested. We defer to 
future work the construction of a new dust map that in- 
corporates the insight gained from these measurements. 

David Schlegel suggested using SDSS stellar spectra to 
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APPENDIX 

EXTINCTION IN DIFFERENT BANDPASSES 

The spectrum-based reddening measurements, in concert with these and the iSchlaflv et al\ (|2010l . SIO) reddening 
measurements, demand that the colors of high lati tude stars with E{B — V) < 1 he dereddened with an Ry = 3.1 



measurements, aemana tnat tne colors or nign lati tuae stars witn — v ) < i De aereddened witn an Uy = d.i 
F99 reddening law rather than an iO'Donnelll ()1994[ ) reddening law. The F99 predictions for Ra-b agree with the blue 
tip values to within 3% in 5 — r, r — i, and i — z, and are off by 6% in u — g. This agreement is sufficiently good that 
we have confidence that reddening in non-SDSS bands can reliably be predicted according to the F99 reddening law. 
Accordingly, we have tabulated the reddening per unit E{B — V^)sfd for a large number of commonly used bandpasses 
(Table [6|). 

For this computation, following the notation of SFD, we compute 



Ami, = —2.5 log 



/ dAW^b(A)5(A)10-'4(^)'^"i^»/2-5 



/ dXWb{X)S{X) 



(Al) 



where W is the system throughput for the band 6, the source spectrum is S in photons/ A/s, A is the extinction law, 
normalized so Ai = 1, and Ami /xm is N times the extinction at 1 /zm according to SFD, following the Ry =3.1 
O'Donnell extinction law assumed by SFD. The factor N here serves as a normalization factor; we use N = 0.78 
following the results of SIO. 

The extinction law is parameterized by Ry = Ay /E{B — V). We have tabulated Ami,/E{B — V)sfo for Ry = 2.1, 
3.1, 4. 1, and 5.1, in the lim it that E{B — t^)sFD is small. The traditional value of Ry in the diffuse interstellar medium 
is 3.1 (jCardelli et a/.|[T989() . which is also what we observe in this work. Values of Ry other than 3.1 are included in 
the table only for convenience. 

For the source spectrum, we use for S{X) a synthetic spectrum from iMunari et al\ (pOOl . with Teff = 7000K, 
\ogZ — —1, and \ogg = 4.5. 

We tabulate values of Amb/E{B — V^)sfd for the filters considered by SFD, using the F99 reddening law. For the 
Landoh UBVRkcIkc, CTIO UBVRkcIkc, Stromgren uhv^y, Gunn griz, Spinrad i?s, UKIRT JHKL', and HST 
WFPC2 bandpasses we follow the procedure detailed in the SFD appendix, combining the filter response with the 
optical and detector throughput as well as the atmospheric transparency at the usual site. The SFD appendix used 
estimated filters for the SDSS from 1994; we now use updated curves taken in June 2001. 

The values of Amb/E{B — F)sfd incorporate a few other changes relative to those in SFD. The CTIO atmosphere 
model SFD used cut off at 3200A, slightly affecting the predicted Amb/E{B — V)sfd in blue bands using that model. 
The source spectrum was also spuriously multiplied by a factor of A, for Wb for which SFD did not use a separate 
detector efficiency curve (UKIRT, Gunn, SDSS, WFPC2, and DSS filter systems). These shortcomings have been 
corrected. 

We additionally tabulate Amb/E{B — V)sf'd for new filter sets that have come to prominence since 1998. We 
have added values of Amb/E{B — y)sFD for the PSl filter complement, using the throughput measurements from 
IStubbs et al\ (|2010D . and for the LSST target fi lter complement, using the "baseline" expected throughput, downloaded 
from the LSST public SVN repository ()Tvsonl 2002). Finally, we have computed Amb/E{B — y)sFD for the wide pass 
filters available to the Hubble Space Telescope's WFC3 and ACS detectors, using the STSDAS package. STSDAS is 
a product of the Space Telescope Science Institute, which is operated by AURA for NASA. 

The IMunari et al\ (|2005f ) synthetic spectra cover the wavelength range 25OOA - IO5OOA. This means that values 
of Amb/E{B — V)sFD for bandpasses that cover regions outside this range use extrapolated source spectra. We 
extrapolate into the infrared following a blackbody spectrum (iS'(A) oc A~'^), which is more blue than the typical 
spectrum and so biases the Amb/E{B - V)sfb somewhat high. This affects the UKIRT JHKL' and HST F105W, 
FllOW, F125W, F140W, and F160W filters. The effect is smaU, however: compared to an S'(A) oc A~^ spectrum, 
Amb/E{B - V)sFD is about 1% larger for the UKIRT filters. 

The CTIO atmospheric transparency spectrum used for the CTIO and Landolt filters cuts off in the blue at 3200A. 
Blueward of this cutoff we extrapolate that the transmission is constant, rendering the transmitted spectrum too blue 
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Table 6 

F99 Reddening in Different Bandpasses 



Bandpass 




2 1 


1 


4 1 


^ 1 


Bandpass 




2 1 


0. 1 


4 1 


^ 1 


Landolt U 


3508.2 


5.614 


4.334 


3.773 


3.460 


WFPC2 F300W 


3087.6 


6.777 


4.902 


4.127 


3.710 


Landolt B 


4329.0 


4.355 


3.626 


3.290 


3.096 


WFPC2 F450W 


4587.0 


4.014 


3.410 


3.132 


2.971 


Landolt V 


5421.7 


2.953 


2.742 


2.645 


2.589 


WFPC2 F555W 


5439.4 


2.976 


2.755 


2.653 


2.594 


Landolt R 


6427.8 


2.124 


2.169 


2.189 


2.201 


WFPC2 F606W 


5984.8 


2.469 


2.415 


2.389 


2.375 


Landolt / 


8048.4 


1.410 


1.505 


1.548 


1.573 


WFPC2 F702W 


6887.9 


1.850 


1.948 


1.994 


2.020 


CTIO U 


3733.9 


5.170 


4.107 


3.628 


3.355 


WFPC2 F814W 


7940.0 


1.452 


1.549 


1.594 


1.620 


CTIO B 


4308.9 


4.382 


3.641 


3.300 


3.104 


WFC3 F105W 


10438.9 


0.981 


0.969 


0.964 


0.961 


CTIO V 


5516.6 


2.857 


2.682 


2.600 


2.553 


WFC3 FllOW 


11169.7 


0.907 


0.881 


0.870 


0.863 


CTIO R 


6520.2 


2.055 


2.119 


2.149 


2.166 


WFC3 F125W 


12335.5 


0.778 


0.726 


0.701 


0.687 


CTIO / 


8006.9 


1.420 


1.516 


1.561 


1.587 


WFC3 F140W 


13692.3 


0.672 


0.613 


0.586 


0.570 


UKIRT J 


12482.9 


0.764 


0.709 


0.684 


0.669 


WFC3 F160W 


15258.3 


0.570 


0.512 


0.485 


0.469 


UKIRT H 


16588.4 


0.502 


0.449 


0.425 


0.411 


WFC3 F200LP 


5515.2 


3.457 


2.958 


2.743 


2.625 


UKIRT K 


21897.7 


0.331 


0.302 


0.288 


0.280 


WFC3 F218W 


2248.3 


12.405 


7.760 


5.956 


5.027 


UKIRT L' 


37772.5 


0.159 


0.153 


0.150 


0.148 


WFC3 F225W 


2394.0 


10.907 


6.989 


5.458 


4.666 


Gunn g 


5200.0 


3.225 


2.914 


2.770 


2.687 


WFC3 F275W 


2742.5 


7.986 


5.487 


4.488 


3.963 


Gunn r 


6628.5 


1.959 


2.055 


2.099 


2.125 


WFC3 F300X 


2934.5 


7.437 


5.228 


4.331 


3.854 


Gunn i 


7898.6 


1.454 


1.555 


1.601 


1.628 


WFC3 F336W 


3366.4 


5.835 


4.453 


3.853 


3.519 


Gunn z 


9050.1 


1.188 


1.234 


1.255 


1.267 


WFC3 F350LP 


5877.1 


2.876 


2.624 


2.509 


2.443 


Spinrad Rg 


6927.3 


1.810 


1.921 


1.972 


2.002 


WFC3 F390W 


3994.8 


4.803 


3.896 


3.481 


3.244 


Stromgren u 


3510.0 


5.539 


4.305 


3.759 


3.452 


WFC3 F438W 


4335.3 


4.347 


3.623 


3.288 


3.095 


Stromgren b 


4670.5 


3.916 


3.350 


3.089 


2.938 


WFC3 F475W 


4785.0 


3.755 


3.248 


3.013 


2.878 


Stromgren v 


4119.2 


4.619 


3.793 


3.411 


3.191 


WFC3 F475X 


4969.7 


3.548 


3.116 


2.917 


2.803 


Stromgren /3 


4861.3 


3.655 


3.183 


2.966 


2.840 


WFC3 F555W 


5302.8 


3.135 


2.855 


2.726 


2.652 


Stromgren y 


5478.9 


2.862 


2.686 


2.605 


2.557 


WFC3 F600LP 


7362.3 


1.688 


1.781 


1.824 


1.849 


SDSS u 


3586.8 


5.419 


4.239 


3.715 


3.419 


WFC3 F606W 


5868.5 


2.581 


2.488 


2.445 


2.421 


SDSS g 


4716.7 


3.843 


3.303 


3.054 


2.910 


WFC3 F625W 


6225.8 


2.230 


2.259 


2.273 


2.281 


SDSS r 


6165.1 


2.255 


2.285 


2.300 


2.308 


WFC3 F775W 


7630.9 


1.533 


1.643 


1.694 


1.724 


SDSS i 


7475.9 


1.583 


1.698 


1.751 


1.782 


WFC3 F814W 


7983.1 


1.441 


1.536 


1.580 


1.605 


SDSS z 


8922.9 


1.211 


1.263 


1.286 


1.300 


WFC3 F850LP 


9149.7 


1.168 


1.208 


1.226 


1.237 


DSS-II g 


4620.6 


3.970 


3.381 


3.110 


2.954 


ACS clear 


6211.1 


2.612 


2.436 


2.356 


2.309 


DSS-II r 


6545.5 


1.991 


2.088 


2.133 


2.159 


ACS F435W 


4348.3 


4.330 


3.610 


3.278 


3.087 


DSS-II i 


8111.0 


1.396 


1.487 


1.530 


1.554 


ACS F475W 


4760.3 


3.787 


3.268 


3.028 


2.890 


PSl g 


4876.7 


3.634 


3.172 


2.958 


2.835 


ACS F550M 


5581.0 


2.754 


2.620 


2.558 


2.522 


PSl r 


6200.1 


2.241 


2.271 


2.284 


2.292 


ACS F555W 


5361.3 


3.031 


2.792 


2.682 


2.618 


PSl i 


7520.8 


1.568 


1.682 


1.734 


1.765 


ACS F606W 


5901.0 


2.555 


2.471 


2.431 


2.409 


PSl z 


8665.3 


1.258 


1.322 


1.352 


1.369 


ACS F625W 


6298.1 


2.171 


2.219 


2.241 


2.254 


PSl y 


9706.3 


1.074 


1.087 


1.094 


1.097 


ACS F775W 


7673.5 


1.520 


1.629 


1.679 


1.708 


PSl ii) 


6240.8 


2.425 


2.341 


2.302 


2.280 


ACS F814W 


8012.4 


1.432 


1.526 


1.569 


1.594 


LSST u 


3693.2 


5.243 


4.145 


3.652 


3.373 


ACS F850LP 


9007.5 


1.196 


1.243 


1.265 


1.277 


LSST s 


4797.3 


3.739 


3.237 


3.006 


2.872 


DES g 


4796.6 


3.739 


3.237 


3.006 


2.872 


LSST r 


6195.8 


2.245 


2.273 


2.286 


2.294 


DBS r 


6382.6 


2.113 


2.176 


2.205 


2.221 


LSST i 


7515.3 


1.571 


1.684 


1.737 


1.767 


DES i 


7769.0 


1.490 


1.595 


1.644 


1.672 


LSST z 


8664.4 


1.259 


1.323 


1.353 


1.370 


DES z 


9108.2 


1.175 


1.217 


1.236 


1.247 


LSST y 


9710.3 


1.075 


1.088 


1.094 


1.098 


DES Y 


9850.4 


1.051 


1.058 


1.061 


1.063 



Note. — Ai,/E{B — V)sfd in different bandpasses b, evaluated according to an F99 reddening law with normalization = 0.78 and 
Ry = 2.1, 3.1, 4.1, and 5.1, using a 7000 K source spectrum. The column A^fj gives the throughput-weighted mean wavelength in the 
bandpass. W hen used with Ry = 3.1, these give the coefficients to use with E(B — V)sfd to get reddenings consistent with the results of 
this work and lSchlaflv et oil II2010I ). The values for other Ry are provided only for convenience. 



and deriving Amii/E{B — F)sfd somewhat too large. This affects only the CTIO and Landolt U bands. 
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